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Abstract 

The chiral extrapolation of the vector meson mass calculated in partially-quenched lattice sim- 
ulations is investigated. The leading one- loop corrections to the vector meson mass are derived 
for partially-quenched QCD. A large sample of lattice results from the CP-PACS Collaboration 
is analysed, with explicit corrections for finite lattice spacing artifacts. To incorporate the effect 
of the opening decay channel as the chiral limit is approached, the extrapolation is studied us- 
ing a necessary phenomenological extension of chiral effective field theory. This chiral analysis 
also provides a quantitative estimate of the leading finite volume corrections. It is found that the 
discretisation, finite-volume and partial quenching effects can all be very well described in this 
framework, producing an extrapolated value of Mp in excellent agreement with experiment. This 
procedure is also compared with extrapolations based on polynomial forms, where the results are 
much less enlightening. 
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I. INTRODUCTION 



There has been great progress in lattice QCD in recent years, associated both with Moore's 
Law and with improved algorithms, which mean that one can work with larger lattice spac- 
ings and still approximate the continuum limit well. The CP-PACS group has devoted 
considerable effort to the study of the masses of the lowest mass baryons and vector mesons. 
This has led, for example, to a comprehensive set of data for the mass of the p-meson in 



partially quenched QCD, with exceptionally small statistical errors We shall exploit this 
data. 

The remaining barrier to direct comparison with experimental data is the fact that cal- 
culations take much longer as the quark mass approaches the chiral limit. Indeed the time 
for a given calculation scales somewhere in the range to m~^, depending on how hard 
one works to preserve chiral symmetry. As a result there has been considerable interest in 
using chiral perturbation theory (^PT), an effective field theory (EFT) built on the sym- 
metries of QCD, to provide a functional form for hadron properties as a function of quark 
mass . In principle, such a functional form can then be used to extrapolate from the 
large pion masses where lattice data exists to the physical value. Unfortunately, there is 
considerable evidence that the convergence of dimensionally regularised xPT is too slow for 
this expansion to be reliable at present 

On the other hand, it can be shown that a reformulation of xPT using finite-range regu- 
larisation (FRR) effectively resums the chiral expansion, leaving a residual series with much 
better convergence properties The FRR expansion is mathematically equivalent to 

dimensionally regularised xPT to the finite order one is working . Systematic errors 

associated with the functional form of the regulator are at the fraction of a percent level 
j^. A formal description of the formulation of baryon yPT using a momentum cutoff (or 
FRR) has recently been considered by Djukanovic et al. [12j. The price of such an approach 
is a residual dependence on the regulator mass, which governs the manner in which the 
loop integrals vanish as the pion mass grows large. However, if it can be demonstrated that 
reasonable variation of this mass does not significantly change the extrapolated values of 
physical properties, one has made progress. This seems to be the case for the nucleon mass 
^ and magnetic moments for example, where "reasonable variation" is taken to be 
±20% around the best fit value of the optimal regulator mass. 

In order to test whether the problem is indeed solved in this way one needs a large 
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body of accurate data. This is in fact available for the p meson, where CP-PACS has 
carried out lattice simulations of partially quenched QCD (pQQCD) with a wide range 
of sea and valence masses. This sector requires a modified effective field theory, namely 
partially-quenched chiral perturbation theory (PQxPT) Formal developments in 

this field have made significant progress in the study of a range of hadronic observables — 



see Refs. 



3ll2 
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19 



21 



{ for example. 

In considering the p meson, analysis of modern lattice results requires one to extend 
beyond the low-energy EFT. Near the chiral limit, the p decays to two energetic pions, 
whereas at the quark-masses simulated on the lattice the p is stable. The pions contributing 
to the imaginary part of the p mass cannot be considered soft, and therefore cannot be 
systematically incorporated into a low-energy counting scheme 0, 24 1. Because almost 
all the lattice simulation points in this analysis lie in the region m^r > 'n^p/2, it is evident 
that the extrapolation to the chiral regime will encounter a threshold effect where the decay 
channel opens. To incorporate this physical threshold, we model the p —>■ tttt self-energy 
diagram constrained to reproduce the observed width at the physical pion mass. Including 
this contribution also provides a model of the finite volume corrections arising from the 
infrared component of the loop integral. In particular, we can also describe the lattice 
results in the region m.,^ < mp/2, where the decay channel is still energetically forbidden 
because of momentum discretisation. 

This large body of pQQCD simulation data is analysed within a framework which incor- 
porates the leading low-energy behaviour of partially-quenched EFT, together with a model 
for describing the decay channel of the p meson. Finite-range regularisation is implemented 
to evaluate loop integrals, for reasons discussed above. The aim is to test whether it pro- 
duces a more satisfactory description of the complete data set than the more commonly 
used, naive extrapolation formulas. A condensed version of some the work featured here has 
been reported in Ref. 

The next section summarises the finite-range regularised forms for the self-energy of the 
p meson in the case of pQQCD. Section IIIII discusses the data used from the CP-PACS 
Collaboration We then give details of the chiral fits in Sec. IIVI Finally, Sec. IVl reports 
the experimental determination of the p meson mass at the physical point. 
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II. SELF-ENERGIES FOR THE PARTIALLY QUENCHED ANALYSIS 



Theoretical calculations of dynamical-fermion QCD provide an opportunity to explore 
the properties of QCD in an expansive manner. The idea is that the sea-quark masses 
(considered in creating the gauge fields of the QCD vacuum) and valence quark masses 
(associated with operators acting on the QCD vacuum) need not match. Such simulation 
results are commonly referred to as partially quenched calculations. Unlike quenched QCD, 
which connects to full QCD only in the heavy quark limit, partially quenched QCD is not 
an approximation. The chiral coefficients of terms in the chiral expansion (such as the axial 
couplings of the vr and 77') are the same as in full QCD. Hence, the results of partially 
quenched QCD provide a theoretical extension of QCD [iS]. QCD, as reahzed in nature, is 
recovered in the limit where the valence and sea masses match. 

In this section we explain the form of the finite-range regularised chiral extrapolation 
formula in the case of partially quenched QCD (pQQCD) — i.e., the case where the valence 
and sea quarks are not necessarily mass degenerate. This work extends on the early work of 
Ref . 3| and the more recent analysis of Ref . , which studied the case of physical (full) 
QCD. This discussion includes the self-energies S^^ and S^^^ (corresponding to Eqs.(3) and 
(4) in Ref. [2^) and includes the self-energy contributions associated with the double hairpin 
diagrams surviving to some extent in pQQCD. We restrict our attention to the case where 
the two valence quarks in the vector meson are degenerate. 

We introduce the following notation. Mps{v){Py ^sea; /^vai' '^vai) I'efers to the pseudoscalar 
(vector) meson mass, with the first two arguments referring to the gauge coupling and sea 
quark mass, while the last two refer to the valence quark masses. Throughout the paper it 
will be convenient to use a short hand notation: 



^non-deg ^ M(/?, /t^ea; /^sea, /^val) 
Mdeg = M(/?, Ksea; «:val, «;val) 
j^umt _ jVf /^gg^- /^gg^^ ^gg^)_ 

where the superscript unit refers to the unitary data with kI^^ = k^^i = Usea] deg refers to 
the "degenerate" data, where n\g^\ = k^^i and these are not necessarily equal to Kgca] non-deg 
refers to the non-degenerate case where k^^^j = Ksea 7^ ^^vai- 

Derivation of the pQQCD chiral expansion for mesons can be described via the diagram- 
matic method Q|, where the role of sea-quark loops in the creation of pseudoscalar meson 
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FIG. 1: The leading nonanalytic contribution to the chiral expansion of the /3-meson mass. The 
meson dressing (a) and its associated quark-flow diagram (b) are illustrated. 




FIG. 2: Two-pion contributions to the positively charged p-meson self-energy. The quark-flow dia- 
grams corresponding to the meson dressings of (a) and (c) are illustrated in (b) and (d) respectively. 



dressings of the vector meson is easily observed. Consider the simplest case of the n-uj dress- 
ing of the p meson depicted in Fig. Q which gives rise to the leading nonanalytic (LNA) 
contribution to the chiral expansion of the p-meson mass in full QCD. Here the positive 
charge state of the p is selected to simplify the derivation. 

The two-pion contribution to the p-meson self-energy is depicted in Fig. |21 This channel 
gives rise to the next-to-leading nonanalytic (NLNA) contribution to the chiral expansion of 
the p-meson mass. Importantly, this contribution also ensures that the p develops a finite 
width as the two-vr decay becomes accessible. 



As these channels can only come about through the inc 



usion of a sea-quark loop, the 



expressions for the pionic self-energies are as given in Ref. 27|, but with the pion mass being 
j^noji-deg^ Corresponding to one valence quark and one sea quark. 
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FIG. 3: The r/' contribution to the p-meson self energy and its associated quark flow diagrams in 
pQQCD. While diagram (c) appears in quenched QCD, it is complemented by an infinite series of 
terms, the first few of which are depicted in diagrams (d) through (f). 



In the case of the partially-quenched rj' contributions to the chiral expansion of pQQCD, 

28^ ) that the Witten- Veneziano coupling 



we simplify the calculation by assuming (as in Ref. 
is simply a constant, — /Uq, and take the limit fio ^ oo after summing all relevant sea quark 
bubble diagrams. The philosophy is that the physical rj' mass is so large that we can ignore 
rj' loops in full QCD. 

The 1]' contributions to the p-meson self energy are depicted in Fig. Efa) with the asso- 
ciated quark flow diagrams illustrated in Fig. Efb) through (f). While Fig. Efc) appears in 
quenched QCD, it is complemented by an infinite series of terms, the first few of which are 
depicted in Fig. IHfd) through (f). Only the sum of Fig. Efb) through (f) and beyond raise 
the 7]' mass in full QCD. 

In partially quenched QCD it is essential to track the masses of the pseudoscalar mesons 
contributing to the quark-flow diagrams of Fig. 01 Whereas Fig. Efb) involves a non- 
degenerate pseudoscalar meson. Fig. Efc) involves degenerate pseudoscalar mesons, and the 
remaining figures involve both degenerate and non-degenerate pseudoscalar mesons. The 
sum of diagrams (b) through (f ) of Fig. IHl and beyond generate an rj' propagator propor- 
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tional to 



92 glf^l 



1- . ..o + 



P5 ^ ^"P5 



Upon summing the terms in [■■■], Eq. takes the form 



2 2 _|_ Afunit'^ 2 



We note that for equal valence- and sea-quark masses where Mp^ = Mp°^~'^'^^ = Mp^^^, 
Eq. (j2I) describes the propagation of a heavy rj' meson. Upon taking fi^ oo 



,tot _ 92 



^tot ^2 V /_ /gN 

By subtracting and adding + Mp'^ j to the first and second terms respectively, this 
can be rearranged as 

(7 , ~ 5 5 1 15 (4) 

(g2 + M^""-^'^^ ) (g2 + M^f ) (g2 + M^f )2 
which clearly vanishes when the sea and valence quark masses are equal. 

The quenched coupling of a light rj' meson (degenerate with the pion) depicted in Fig. 
is complemented by an infinite sum of sea-quark loop contributions in pQQCD, the first few 
of which are depicted in Fig. IH^b) through (d). Using similar techniques, one can show that 
the sum of these graphs leads to the propagation of a heavy r]' and upon taking fio oo the 
contribution of Fig. ^ vanishes. 

We note that the coupling g2 introduced in Ref. Q], takes the value 0.75, where in 
terms of the usual upn coupling constant, g^jpn = 16 GeV~^, (7I = g'^p^f^/i. Introducing 
fpTTu: ~ t^p92 (with yUp(7r) the physical pM mass), the total, non-trivial rj' contribution can be 
written in the form (following Ref. [27[ for large vector meson mass): 



y^P _ f p-TTW f k U jk) dk // ,^non-(ieg\2 _ ( T\/rdeg\2\ 



In summary, we find 

T.TOT = + + K'P ' (6) 
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FIG. 4: The quenched couphng of a hght r]' meson (degenerate with the pion) depicted in (a) is 
complemented by an infinite sum of sea-quark loop contributions in pQQCD, the first few of which 
are depicted in (b) through (d). 



where the pion self-energies are given by: 



yp _ ffmu) f°° k^ "^Iwik) dk 

" 3n^f^l uj^{k)Mk) + ^M^p)' 



and 



coKk) = + (m;°"-'^^)2 . (9) 

As the LO meson contributes via a sea-quark loop as illustrated in Fig. ^ the mass splitting 
between the p and u is 

AM^p = M^'^-'^'^ - Mi'^ . (10) 

We note that {ujT^{k) + AM^p) > for all quark masses and nontrivial momenta considered 
in the lattice analysis. We use the values /ptttt = 6.028 (obtained from Ref. I^l) and = 
92.4 MeV. We use a standard dipole form factor, which takes the form 

= (A2 + p)2 ' 

u^i^{k) = u{k), 



where the second factor in m^^(/c) ensures the correct on-shell normalization condition. 



To account for finite volume artefacts, the self-energy equations are discretised so that 



only those momenta allowed on the lattice appear |27|, |3 
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k^dk 



d^k 



1 
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k 



2vr(i, j, k) 



x,y,z 



'ID 



(12) 



x,y,z 



The purpose of the finite-range regulator is to regularise the theory as k^, ky, k^ tend to 
infinity. Indeed, once any one of kx-, ky or k^ is greater than ~ lOA the contribution to the 
integral is negligible and thereby ensuring convergence of the summation. Hence, we would 
like the highest momentum in each direction to be just over lOA. For practical calculation, 
we therefore use the following to calculate the maxima and minima for i, j, k: 

'lOA a 



27r 
lOA a 



{x,y,z) 



27r 



N[x,y,z) 



+ 1, 
-1, 



where [. . .] denotes the integer part. 

We study a range of values of A, starting with the value, A^^^ = 630 MeV, used in Ref. |27|. 
Figure El shows the various self-energy contributions, and S^,^ as a function of 

j^non-deg ^^^^ -g^^^ ^ © and © ) for the representative (/?, K^ea) = (2.10,0.1382) dataset. 
In Sec. II V CI we perform a highly constrained fit to a large "global" dataset, and this enables 
us to determine a best value of A which minimises the global x^. 



III. OVERVIEW OF CP-PACS DATA 

In Ref. Q], the CP-PACS collaboration published meson spectrum data from dynamical 
simulations for mean-field improved Wilson fermions with improved gluons at four different 
/3 values. For each value of /3, ensembles were generated for four different values of Ksea 
- giving a total of 16 independent ensembles. Table E] summarises the lattice parameters 
used. The pseudoscalar- to vector-meson mass ratio Mp™*/^y"** (which gives a measure of 
the mass of the sea quarks used) varies from 0.55 to 0.80 and the lattice spacing, a, from 
around 0.09 to 0.29 fm. For each of the 16 ensembles there are five Ky^i values considered Q|. 
Thus there are a total of 80 {My^,Mp^) data points available for analysis. In our study 
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FIG. 5: The self-energy contributions (see Eqs. EllHl&ISl) versus Mp°g '^'^^ data for the ensemble 
(/3,Ksea) = (2.10,0.1382) using A = 630 MeV. 



we consider the two cases where the lattice spacing is set using either tq [31I] or the string 
tension a. 

In the absence of the full set of original CP-PACS data, we generate 1000 bootstrap 
clusters for all Mps and My data using a Gaussian distribution whose central value and 
full width half maximum (FWHM) are the same as the central values and errors published 
in the Table XXI of Ref. P|. Of course, our errors are totally uncorrelated throughout 
- i.e., each Mv(/3, Ksea; «^vab ^vai) bootstrap cluster is uncorrelated with the corresponding 
Mps{l3 1 f^sea'i f^la\, f^lai) bootstrap clustcr. Furthermore, the M{(3, Kgea.'-, ^^vab '^vai) data is not 
correlated with the M(/?', k^^^; k^^^j, k^^j) data, and nor is the M(/3, Kgea; /^vai; '^vai) data corre- 
lated with that for M(/3, Ksea; ^^vli' ^vli)- 

We expect therefore that the statistical errors of our final results will be overestimates 
of the true error, since we have not benefited from the partial cancellation of statistical 
errors which occurs when combining correlated data. We can obtain a rough estimate of the 
increase in our errors due to the fact that we don't maintain correlations as follows. The 
ratio M]^^*7^r** listed in Table HI is obtained from our bootstrap data. Comparing this 
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TABLE I: The lattice parameters of the CP-PACS simulation used in this data analysis, taken 
from Ref. ^|. The superscript unit refers to the unitary data (i.e., where k^^^j = k^^^j = Kgea)- Note 
that the errors reported in this table are obtained with our bootstrap ensembles (see text). 

with the Mp™*/^y™* data in Table XXI of Ref. fl| (which benefits from the cancellation of 
correlations), we see that a very rough estimate of the effect of ignoring correlations is to 
increase the errors by ~20%. We expect that a similar increase in errors will apply to other 
quantities. 

The lattice spacings a^^^a were obtained from Table XII of Ref. Q| using rg = 0.49 fm and 
a/o" = 440 MeV. Again we generated 1000 bootstrap clusters with a Gaussian distribution, 
as in the meson mass data. Figure El shows the unitary (i.e., k^^^i = K^g^j = Kgca) pseudoscalar 
mass plotted against the lattice spacing, a^g , for the 16 ensembles in Table |l] (Note that 
(Mp™*)^ is a direct measure of the sea quark mass since, from PCAC, (Mp™*)^ oc m1^^.) 
Also shown for reference are the physical pseudoscalar mesons vr, K and "?7s". Note the large 
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range of both the lattice spacing a^o and ml^^ in the simulations, and that the lattice spacing 
is primarily determined by the (5 value, rather than the value of ml^^. 

The physical volume for these 16 ensembles is La 2.5 fm for the (3 = 1.80, 1.95 and 2.10 
cases, but the (3 = 2.20 ensemble had a slightly smaller physical volume. The associated finite 
volume effects are incorporated through evaluating the chiral loops by explicitly summing the 
discrete pion momenta allowed on the lattice. We treat all 16 ensembles on an equal footing. 
The finite-volume effects are corrected when making contact with the physical observables 
by evaluating the chiral loop integrals with continuous loop momenta. 

The action used in Ref. Q] is mean-field improved, rather than non-perturbatively im- 
proved and will therefore have some residual lattice systematic errors of 0(a) Q- We fit 
the data assuming both 0{a) and O(a^) effects in sections llV Bl and IIV CI and are thus able 
to obtain continuum predictions. Our empirical analysis suggests that nonanalytic terms 
generated in a dual expansion of both a and are either small or 

can be absorbed at present into the 0{a) and O(a^) effects considered here. 



IV. FITTING ANALYSIS 



A. Summary of Analysis Techniques 

Our chiral extrapolation approach is based upon converting all masses into physical units 
prior to any extrapolation being performed. An alternative approach would be to extrapolate 
dimensionless masses (in lattice units) |l!] . However, using physical units offers the following 
key advantages: 

• The data from different ensembles can be combined in a global fit. If the masses were 
left in dimensionless units, there would be no possibility of combining data at different 
lattice spacings together into such a global fit. 

• Dimensionful mass predictions from lattice simulations are effectively mass ratios, and 
hence one would expect some of the systematic (and statistical) errors to cancel. That 
is, M''^""^"' = M* X a~^ = M*/M* x M^^\ where Vt is the quantity used to set the 
lattice spacing, a, the superscripts expt refer to the lattice mass estimate and the 
experimental value, respectively. 

• Lattice data are not in the regime where the coefficients of nonanalytic terms in the 
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FIG. 6: The range of sea quark mass (Mp^**)^ and lattice spacing, a^Q, covered by the CP-PACS 
data displayed in Table (Mp™*)^ is the pseudoscalar meson mass squared at the unitary point; 
i.e., where Kvai = «^sea- The experimental points for the tTjK and "r/^" mesons are also shown for 
reference. 

chiral expansion can be reliably constrained. These dimensionful quantities must there- 
fore be fixed to their phenomenological values. 

In this paper we use two different quantities for setting the scale, tq and a - although we 
have a preference, as explained in Sec. II for tq. Table |I] lists values for a^^ and a^- 

In our chiral extrapolations we use two basic fitting functions, the finite-range regulari- 
sation method (hereafter referred to as the "Adelaide" method) 



where T^j-ot is calculated using Eq. ©, and a naive polynomial fit 



(13) 



(14) 
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In each case we refer to these fits as "cubic" since they include cubic terms in the chiral 
expansion of rriq oc (Mps)"^. We also perform fits with the coefficient Og set to zero in 
Eqs. (jl3p and (jl4j) . referring to these as "quadratic". It is worth noting that the dominant 
functional form of My with (Mpg)^ is linear (see, for example, Fig. [71 where the LHS of 
Eqs. fll3|) and ()14|) are referred to as My^). We exploit this fact in the fitting functions given 



above. In particular, this is why the Adelaide fit uses y {My^Y ~ ^tot on the LHS rather 
than (Mf - 

TiTOT which would, a priori, be an equally valid chiral expansion. Thus, with 
the above functional forms, we expect the coefficients, a^, to be small for n > 4 and this is, 
in fact, what we find. 

In the following, we first use Eqs. (fT^ and (fT^ to fit the 16 ensembles in Tahle^separately. 
We then turn to a global fit, where all 16 ensembles are combined in one fitting function. 

B. Individual ensemble fits 

We begin our analysis by fitting the meson spectrum of each of the 16 ensembles listed 
in Tahle^ separately. In this section we use ro to set the scale and select a value of A = 650 
MeV (see Sec. HVTT|l . We use both the Adelaide, Eq. and naive, Eq. dUl), fitting 

functions and restrict our attention to quadratic (ag = 0) chiral fits, since there are only 
five {My^^ , Mp^) data points available for each ensemble. The coefficients, 00,2,4, obtained 
by fitting My against Mps with both the naive (Eq. (Hlj)) and Adelaide (Eq. (fT^ ) fitting 
functions are listed in Table |nj We see that the 04 coefficients are both small and generally 
poorly determined, confirming our decision to fit to the quadratic, rather than the cubic, 
chiral extrapolation form. We note also that there is some agreement between the naive and 
Adelaide ao,2 coefficients, although their variation with K^en tends to be different. 

Figure [7| shows the results of these fits for the (/5, Usea) = (2.10, 0.1382) ensemble, which is 
a good representative of all of them. The scale is set from tq, which is our preferred method 
(see Sec. IIVC|) . Note that this ensemble's {a,mj^J coordinates are closest to the physical 
point (a,m^gg^) = (0,m„^rf) for those ensembles with La ~ 2.5 fm - see Fig. IHl 

The values of ao.2 in Table ITTl hint at a systematic variation of ao_2 with a^g. To check this 
we plot ao and 02 against (for both the linear and Adelaide fits) in Figs. |H1 and IHl This 
motivates a continuum extrapolation of the form 





-individual 
0,2 



a, 



(15) 



14 



/3 


^sea 


„ naive 
"0 

[GeV] 


„adel 
"0 

[GeV] 


„naive 
"2 

[GeV-1] 


„adel 
"2 

[GeV-1] 


„naive 
0-4 

[GeV-3] 


^adel 

[GeV~3] 


1.80 
1 80 
1.80 
1.80 


0.1409 
1 4S0 
0.1445 
0.1464 


71 2+1"* 

0.73^2 

0.721^ 


0.701^ 

0.756ti^ 
0.769^^3 


0.461^ 
48+s 
0.431^ 
0.491^ 


0.541^ 
51 +5 
0.441« 
0.431^ 


-0.011? 

0.011^ 
-0.0216 


-0.091^ 
-0 08+6 
-O.Oll^ 
0.0071g 


1.95 
1 9.5 
1.95 


0.1375 
1 390 
0.1400 


76+2 
76+2 

U. ( U„2 

0.785ti2 


75+2 
772+1'' 
0.803tii 


49+'* 
0.43tl 


0.53+4 
49+* 
0.4411 


-0.05+3 
-0.0113 


-0.08+3 
-0 05+* 
-0.0211 


1.95 


0.1410 


0.766t« 


0.799^1^ 


0.481^ 


0.4511 


-0.0311 


-0.0311 


2.10 
2.10 
2.10 


0.1357 
0.1367 
0.1374 


829+^'' 
n 794+" 
0.807ti| 


820+1'' 
797+11 

o.822t;;| 


42+5 
0.50+3 
48+f 


46+5 
0.53+3 
49+1 


-0 02+3 
-0.06+3 
-0.0511 


-0.05+1 
-0.08+2 
-0.061^ 


2.10 


0.1382 


0.781tjo 


0.814^10 


0.5311 


0.501^ 


-0.0812 


-0.0712 


2.20 


0.1351 


0.84+3 


84+^ 


43+** 


46+* 


-0.02+^ 


-0 04+6 


2.20 
2.20 


0.1358 
0.1363 


0.83^2 
0.801^ 


0.84^2 

0.811:^ 


0.441^ 
0.511* 


0.461^ 
0.521* 


-0.031^ 
-0.071« 


-0.051^ 
-0.0816 


2.20 


0.1368 


0.781^ 


0.8012 


0.5211 


0.511^ 


-0.0616 


-0.0616 



TABLE II: The coefficients obtained from fitting My data against Mp^. using both the naive and 
Adelaide fits (i.e., Eqs. ((TH) and (fT^ . respectively) for each of the 16 ensembles listed in Tabled 
As discussed in the text we restrict these fits to quadratic rather than cubic chiral functions (i.e., 
ag = 0). The scale was set using tq. 

The results of these fits are displayed in Table IIIII From the values of Xq^2™^'^^"''' 
Table ITTTl we can confirm statistically significant 0{a) effects in but not in 02- This will 
be important in determining our fit procedure in Sec. IIV CI 

C. Global fits 

We now turn our attention to an analysis of the whole (degenerate) data set. Figure fTUl 
shows all of the degenerate CP-PACS data with the physical scale set using vq. Since there 
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Mp/ [GeV'] 



FIG. 7: A plot of M^"^ versus Mps data for the ensemble (/5,Ksea) = (2.10,0.1382) together 
with the results of the quadratic Adelaide, Eq. H13|) . and naive, Eq. ((Ill), fits. M^"'' is defined 



as M^"^ 



{My^Y ~ "^TOT for the Adelaide fit (i.e., the LHS of Eq. ^J- note that Stot is 
negative), and M^"* = My^ for the naive fit (i.e., the LHS of Eq. 





^.cont. 

"o 

[GeV] 


vindividual 
^0 

[GeV/fm] 


xl/d.o.f. 


„cont. 

[GeV-1] 


^individual 

[GeV-Vfm] 


xl/d.o.f. 


Naive- fit 
Adelaide- fit 


0.8611^1 
0.8731i° 


-0.531^ 
-0.511^ 


21 / 14 
16 / 14 


0.5111 
0.501^ 


-0.21^23 
-0.06111 


8/14 
10 / 14 



TABLE III: The coefficients obtained from the continuum extrapolation of both the naive and 
Adelaide fits, using the values of ao,2 given in Table Hll with Eq. (fT3)). 

are 16 ensembles with five [My^, ^Ps) values in each (see Sec. IIIIIl . this global fit contains 
80 data points. We expect that this large number of data points should produce a highly 
constrained set of fit parameters, ao,2,...- 

Referring to the coefficients listed in Table |n] and the discussion in the previous section, 
we observe a significant variation amongst the Oq values with lattice spacing, whereas the 
values of 02 are approximately constant with lattice spacing. We also recall that the 04 
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coefficient was undetermined. This suggests that we should allow for some variation of the 
qq coefficient with lattice spacing and in consequence we adopt the following modified version 
of the Adelaide and naive ffiting functions, based on Eqs. (fT^ and (fT^ . 

^[{M^^Y^T^T = (ar* + X^a + X^a') + a^iM^'f + a,{MplY + a,{Mpl)\ (16) 

and 

Mi'' = (ar* + X,a + X^a^) + a^iMplf + a^{Mplf + a^{Mplf , (17) 

respectively. As in the previous section we refer to the above fits as "cubic", since they 
include the ag term (which is proportional to m^). As above we also perform fits with ag 
set to zero, referring to these as "quadratic". These ffiting functions, Eqs. (fT^ and (fTTj) . 
have both 0{a) and O(a^) terms, because the lattice action used is only mean-field improved 
and will contain O(a^) errors, together with some residual errors at 0{a). In the following 
we will experiment by turning off the 0{a) term (i.e. by setting Xi = 0) in Eqs. (fT^ 
and p7|) in order to see whether these residual 0{a) errors are significant. Note that we 
also included 0{a, o?) terms in 02 (and even 04) as a check but found that these ffis were 
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FIG. 9: A continuum extrapolation of the a2 coefficient obtained from both the Adelaide and naive 
fits using Eq. (fT3)) . 



unstable, confirming the findings of the previous section that there are discernible lattice 
spacing effects only in the oq coefficient. 

The global fits used both tq and the string tension, a, to set the scale. Thus we have a 
large number of fit types which are summarised in Table IIVI Indeed, there are two choices 
from each of the four columns in Table IIVI making a total of 2^ fitting procedures. In the 
following all 2^ fits were performed. As noted above, the global fits contain 80 data points, 
and the largest number of fitting parameters is six (a™"*, Xi, X2, ^4 and a6)j3|. Thus 
the global fits are highly constrained. 

Before presenting results from the global fits, we recall our discussion of the A parameter 
in Sec. |n] The Adelaide approach motivates the introduction of the mass scale. A, as 
corresponding to the physical size of the pion source in the hadron which controls the chiral 
physics. It serves to separate the region where chiral physics is important from that where 
the internal structure of the hadron, which is not part of the effective field theory, becomes 
dominant. Since we are performing a highly constrained fit procedure, we are able to derive 
the best A value from the data as follows. Figure ITTl shows the x^/d-o.f. as a function of A 
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FIG. 10: Illustration of the 80 lattice data points of the degenerate CP-PACS data set, with the 
physical scale of a^g. 



Approach 


Chiral Extrapolation 


Treatment of Lattice 


Lattice Spacing 






Spacing Artefacts in oq 


set from 


Adelaide 


Cubic 


ao term has 


ro 


i.e. ealT^ 


i.e. 0{M^g) included 


0{a + a^) corrections 




Naive 


Quadratic 


ao term has 


a 


i.e. ealTTl 


i.e. no 0{M^g) term 


only 0(0^) corrections 





TABLE IV: The different fit types used in the global analysis. Fits for each of the 2^ choices 
depicted above were performed. 

and we see that for the results where rg was used to set the scale, all four fitting types display 
near identical behaviour, with a distinct minimum at A ^ 650 MeV. In other words, the 
behaviour is independent of whether we chirally expand to 0{Mpg) or 0{Mpg) or whether 
we allow for lattice systematics in the ao coefficient of 0{a + a^) or 0{a^). 
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For the case where a is used to set the scale, Fig. shows that there is a distinct 
minimum at A ^ 550 MeV for all four fitting types. All other things being equal, there is no 
preference between the cubic or quadratic chiral fits. However, the plot clearly shows that 
the 0{a + a?) fits are favoured, in that they produce a lower value than the pure 0{a^) 
fits. As a test, we have also fitted the data with a pure 0{a) correction (i.e., Eq. (fTBj) with 
X2 = 0) and found that the values for this fit overlay those from the 0{a + a'^) fits. This 
is strongly suggestive that the dominant lattice-spacing systematic is 0{a) when the string 
tension is used to set the scale. 




400 500 600 700 

A [MeV] 



FIG. 11: Variation of the x'^/d-o.f against A. The insert shows a closeup of the minimum for the 
quadratic fit with 0{a'^) errors using ro to set the scale (i.e. the preferred fitting procedure). The 
dashed horizontal line represents an increase of x'^ from its minimum value by unity for this fit 
procedure (i.e. it represents one standard deviation) - see Sec. The intercept of this dashed line 
with the x^ curves (at A =620 and 690 MeV) is used to derive the range of acceptable A values. 

We now use these values of A (550 MeV and 650 MeV for the a and ro cases, respectively) 
to perform the 16 global fits discussed above and listed in Table IIVI Table |V| gives the 
coefficients and the x^ /d.o.f. for each fit. 
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Fit 


Scale 


ag°'^* Xi X2 


02 


a4 


ag 


X^/d.o.f. 


Approach 


from 


[GeV] [GeVfm-i] [GeVfm-^ 


] [GeV-i] 


[GeV-3] 


[GeV-^l 

L J 








Cubic chiral extrapolation; contains 0{a + a?) 






Adelaide 


rn 
' u 


QMA+_Xl -0.111}^ -l.ltl 


0.4711 


-0.021?„ 


-0.0211 


38 / 74 


A d pi a i d p 


(J 


0.836t?i -0.37ti« -0.212 


0.4411 


0.041^ 


-0.06+1 


53 / 74 


Naive 


u 


0.81911? -0.151^5 .1.1+3 


0.561« 


-o.ieifo 


0.0511 


77 / 74 


Naive 




0.805j:i« -0.38111 -0.31^ 


0.5711 


-0.181?o 


0.061^ 


73 / 74 






Cubic chiral extrapolation; ao contains 0{i 


a^) only 






Adelaide 


ro 


0.8351« — -1.4011 


0.4811 


-0.031?o 


-0.0211 


39 / 75 


A d pi a i d p 


(J 


0.8071® — -1-241^ 


0.4311 


0.061^ 


-0 06+'^ 


67 / 75 


Naive 


ro 


0.8061^0 — -1.491^ 


0.561^ 


-0.171?o 


0.0611 


78 / 75 


Naive 


(J 


0.7751g — -I.3II4 


0.5611 


-0.161?o 


0.051^ 


87 / 75 






Quadratic chiral extrapolation; 


ao contains 0{a + a?') 






Adelaide 


rg 


0.8401i'^ -0.1111* -1.11^ 


0.4931i? 


-0.0611« 




38 / 75 


A d pi fi i d p 


0- 


0.8291« -0.3711" -0.21^ 


0.4901i? 


-0.05211? 




54 / 75 


Naive 


rg 


0.828111 -0.1611^ -1.11^ 


0.5051i? 


-0.0681?o 




78 / 75 


Naive 




0.8121« -0.37111 -0.31^ 


0.5231i^ 


-0.075111 




74 / 75 






Quadratic chiral extrapolation; 


ao contains 0{a?) only 






Adelaide 


rn 

u 


0.8321^ — -1.4011 


0.4941i? 


-0.0611* 




39 / 76 


Adelaide 


a 


0.79911 — -1.231^ 


0.4861i? 


-0.04611? 




68 / 76 


Naive 


rQ 


0.8151* — -1.491* 


0.5061i? 


-0.0681«o 




79 / 76 


Naive 


a 


0.78111 — -1.3111 


0.52011^ 


-0.069111 




88 / 76 



TABLE V: The results of the global fit analysis. Fits for all 2^ combinations depicted in Table HVl 
are shown. 

Summarising the results of Table IVl and Fig. ^2 (and referring to the 2^ different fit types 
listed in Table HV|) we note: 

Fit Approach: The Adelaide method always gives a smaller than the naive approach, 
confirming it as the preferred chiral extrapolation procedure. 

Chiral Extrapolation: The cubic chiral extrapolation (i.e., introducing the 0{Mpg) term 
in Eqs. (|T?)|l and ()17|1 ) leads to a poorly determined ag coefficient in all cases. Furthermore it 
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causes the 04 coefficient to become much more poorly determined than it is in the quadratic 
chiral extrapolation cases. 

Treatment of Lattice Spacing Artefacts in From Table |3 when the scale is set from 
ro, the ao and 02 coefficients (which are the dominant terms in the chiral extrapolation) do 
not depend on whether (9(a + a^) or (9(a^)-only corrections are applied to the coefficient. 
Indeed, when O(a^) only corrections are used, the error in oq is reduced. However, when the 
scale is set from a, ao does depend on how the 0{a) systematics are treated. Note also that 
Xi (i.e. the 0{a) coefficient) in the a ffis are 2-3 times larger than those from the ro ffis. 
This supports our earlier comments above regarding the probable 0{a) systematics when 
the a scale was used. 

Quantity used to set Lattice Spacing: In the Adelaide approach, setting the scale from tq 
gives a significantly smaller than using a (see Fig. Given this, and the comments 
above regarding the probable 0{a) systematics in the a data, we use ro as our preferred 
method for setting the scale. In the naive case the data does not favour setting the scale 
from either ro or cr. 

On the basis of these results, we choose the quadratic chiral extrapolation method with 
the scale set from ro and O(a^) corrections in the ao coefficient to define the central value 
of both the Adelaide and naive fitting procedure. The spread from the other fitting types is 
used to define the error. In Sec. I3 we determine the predictions for physical p— meson mass 
from these fitting types. 

V. PHYSICAL PREDICTIONS 

We are now ready to estimate Mp in the continuum using the Adelaide and naive fits 
performed in the previous section. We obtain this prediction from Eqs. ()13|) and ()14p by 
setting Mpg = Mp°g~'^^^ = Mp™^ = fi^,. We set AM^p to zero in this calculation and also 
note that S^,^ vanishes, as required. Since we are predicting the continuum value for the 
vector meson mass, we calculate the integrals in Eqs. (jZj) and (jH)) directly, rather than using 
the lattice interpretation of the integral in Eq. (fTTjl . 

We obtain continuum estimates of all 2^ fitting types (see Table HV|) using the coefficients, 
^cont a2,4,6, of the global fits of Sec. llVCJl fi.e.. those in Table IV|) . Table EH displays these 
mass predictions. For the Adelaide case we use A = 650(550) MeV, with the physical scale 
set using ro(cr) (see Sec. IIV(J|) . Note that we use the global analysis fSec. IIV C]) . rather than 
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the analysis of Sec. lIVBt which treated the ensembles separately, since the global fit is much 
more tightly constrained. 

It is interesting to study the variation in the prediction of the physical mass of the p with 
A. Figure IT^ shows how the Mp prediction varies with A for each of the Adelaide fits based 
on the ro scale. Using the plot in Fig. we can estimate the range of acceptable A 
values defined by increasing by unity from its minimum, which represents one standard 
deviation. The horizontal dashed line in Fig. ^2 lies along values one more than the 
minimum for the ro case. From this, we determine that the acceptable range of A values is 
620 < A < 690 MeV (see the insert graph in Fig. [TT| . This range of A values is depicted by 
the vertical dashed lines in Fig. ^1 
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FIG. 12: Variation of the physical mass of the p-meson, Mp, as a function of A, using the Adelaide 
approach with the scale set from tq. The two vertical dashed lines define the range of acceptable 
A values (620 MeV < A < 690 MeV) obtained by increasing by unity in Fig. ^2 

From a detailed examination of Table IVll and Fig. we draw the following conclusions: 



The (statistical) errors in the mass estimates are typically around 1%. 
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TABLE VI: Estimates of Mp obtained from the global fits. The Adelaide fits used A = 
650(550) MeV when the scale is set from ro((T). 

• The Adelaide fitting procedure is very stable when the scale is set using tq. When 
the scale is taken from the string tension, the four Adelaide fits are not in mutual 
agreement. The probable reason for this, as outlined in Sec. IIV CI is that using a to 
set the scale introduces significant 0{a) errors. 
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• The Adelaide fitting procedure is quite accurate - lying at most at twice the statistical 
standard error from the experimental value for the tq case. It would require an un- 
certainty of only around 1-2% in tq (and around 2-6% in ^/a) for the Adelaide central 
values to be in agreement with the experimental value of Mp. 

• From Fig. El the variation of Mp with A is small - roughly the same order as the other 
uncertainties. 

• The naive fitting procedure has both a larger spread of values and is further from the 
experimental value than the Adelaide procedure. 

As a result, we conclude that the Adelaide procedure represents a significant improvement 
over the naive approach. 

We obtain final estimates of Mp by taking the the quadratic chiral extrapolation, the 
scale set from r^, and O(a^) corrections in the coefficient in both the Adelaide and naive 
fitting procedure. (The reason for this choice of fit type is described in detail in Sec. IIV CI ) 
The central value in the Adelaide case was obtained at A = 655 MeV from Fig. (which 
is an adjustment of 1 MeV from the value obtained at A = 650 MeV in Table IVHl . The 
value of A = 655 MeV was used since it is where the minimum of occurs in Fig. ^2 We 
obtain an estimate of the error in the fit procedure from the spread in the mass predictions 
using tq for the scale (since we have reservations about the method when the string tension 
is used to set the scale). An estimate of the uncertainty associated with the determination 
of A comes from varying by unity as described above — i.e. by reading off this error from 
the vertical dashed lines in Fig. El Finally then we are led to the following result for the 
physical mass of the p-meson: 

^Adelaide ^ 778 (4) t «+8MeV, (18) 



p 



^Naive ^ 825{A)tfMeV, (19) 



where the first error is statistical, the second is from the fit procedure, and, in the Adelaide 
case, the third error is from the determination of A. The second error on the Adelaide result 
also includes an uncertainty from the choice of finite-range regulator, which contributes 
+3/— 6 MeV to the error j25|. Note that we have not included an error from the determination 
of tq itself. The only other effect which separates our analysis from nature is that the data 
we are analysing contains only 2 rather than 2+1 light dynamical ffavours. We have no way 
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to estimate the residual systematic error from this source once ro is matched to the physical 
value. 

The fit parameters shown in Table |3 allow one to shift each of the simulation results 
to the infinite- volume, continuum limit and to remove the effects of partial quenching — 
hence restoring unitarity in the quark masses. This can be achieved by first considering the 
preferred form of Eq. ()16|) . where Xi = = 0, 



(20) 



This can then be used to let My^^ My"-'^^ at the same time as removing volume and 



discretisation artifacts. Rewriting Eq. ()20|) in terms of the fit vector meson mass gives, 

1/2 



(21) 



In the physical continuum limit, this becomes 



M^™*(a ^ 0) 



.1/2 



• (22) 



With the unknown parameters determined from the best fit to the entire data set, this then 
provides a prescription for restoring the physical limit of the data. Specifically, each of the 
vector meson masses (calculated at finite a and L) are shifted by an amount 

1/2 



SM 



1/2 



(23) 



Here it is understood that the pion masses that one should therefore plot against on the 
X-axis are the unitary pion, at the point where the sea mass is held fixed and the valence 
mass is changed, ie. M(/3, Ksea! /^vah ^^vai) ^ ^(Z^, ^^sea; ^sea, ^^sea)- The final estimate of the 
physical point is given by 



unit I 



0;L 



oo 



estimate 



deg 



+ SM. 



lattice 



(24) 



The results of these shifts are displayed in Fig. El where we observe a remarkable result. 
The tremendous spread of data seen in Fig. is dramatically reduced, with all 80 points 
now lying very accurately on a universal curve. 

One feature of Fig. is the very steady approach to the chiral limit from the physical 
point. This yields a value for the vector meson mass in the chiral hmit (M°) that IS very 
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FIG. 13: The same 80 lattice data points as in Fig. 1101 after correction to restore the infinite- 
volume, continuum and quark-mass unitarity limits. The central curve displays the best-fit from 
the global analysis with A = 655 MeV. The dashed curves show the (statistically constrained) 
bounds on the FRR scale, 0.620 < A < 0.690 GeV. 



near the physical value. To remove the overall scale, we report on the mass shift between 
the physical and chiral values, finding 

Mp-Ml = 3.7(2) (4) (8) MeV, (25) 

where the sources of uncertainty are the same as those reported above for Mp in Eq. (jl8j) . 
This small number indicates that the sigma term for the p is significantly smaller than the 
nucleon, where the derivative is observed to increase near the chiral limit j^. The reduced 
slope in the case of the p arises from the presence of substantial spectral strength in the 
low-energy two-pion channel below the p- meson mass 
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VI. CONCLUSIONS 



In summary, we have tackled the ambitious task of producing a single, unified chiral fit 
to all of the accurate CP-PACS data for the mass of the p-meson in partially quenched 
QCD - i.e., the case where Kgea, 7^ ^vai- As well as using a naive polynomial fit in Mpg, 
we have generalized the Adelaide approach, developing a unified analysis approach, to fit 
the data. This approach enables one to account for finite volume errors by evaluating the 
chiral self-energy contributions with the same momentum discretisation implicit in the lattice 
simulations. In addition, we have been able to quantify the residual 0{a) effects and hence 
carry out a continuum extrapolation (c.f. Figs. |H1 and IHl). 

The obtained for the Adelaide fit, with the physical scale set using ro, is a factor of two 
lower than that for any other method. This provides considerable confidence in the method, 
even before it is used to produce the physical mass of the p. The quality of the fit leads to 
very small (statistical) error bars — see Table IVII 

In addition, it is possible to estimate the systematic errors in the extrapolation to the 
physical p— mass associated with the fitting procedure (both the chiral and continuum- limit 
fitting procedures). In particular, the finite-range regulator parameter. A, is constrained by 
the model-independent lattice QCD data, and the variation of the Mp prediction within this 
range was found to be 1%. 

The curve through Fig.^ldisplays the determined variation of the p- meson mass with pion 
mass. This curve also presents an extrapolation to the physical point, allowing extraction of 
the physical p-meson mass 

Mp = 778(4)li«l«MeV, (26) 

where the first error is statistical, the second is from variations of the fit procedure and the 
third from the determination of A. Whereas the naive fitting procedure leads to a value that 
is 50-60 MeV too high, the result from the chiral analysis is in excellent agreement with the 
experimentally observed mass. 

Acknowledgements 

CRA and WA would like to thank the CSSM for their support and kind hospitality. WA 
would like to thank PPARC for travel support. The authors would like to thank Stewart 
Wright and Graham Shore for helpful comments. This work was supported by the Australian 



28 



Research Council and by DOE contract DE-AC05-84ER40150, under which SURA operates 
Jefferson Laboratory. 



A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 65, 054505 (2002) [Erratum-ibid. 
D 67, 059901 (2003)] arXiv:hep-lat /Ol OSOlSl . 

D. B. Leinweber, A. W. Thomas and R. D. Young, Phys. Rev. Lett. 92 (2004) 242002 
.arXiv:hep-lat /0302020 ' . 

M. Procura, T. R. Hemmert and W. Weise, Phys. Rev. D 69, 034505 (2004) 
|arXiv:hep-lat/0309020| . 

D. B. Leinweber, A. W. Thomas, K. Tsushima and S. V. Wright, Phys. Rev. D 61, 074502 
(2000) arXiv:hep-lat /9906027| . 



S. Durr, Eur. Phys. J. C 29, 383 (2003) arXiv:hep-lat/020805ll 



C. Bernard, S. Hashimoto, D. B. Leinweber, P. Lepage, E. Pallante, S. R. Sharpe and H. Wittig, 
Nucl. Phys. Proc. Suppl. 119, 170 (2003) arXiv:hep-lat/02Q9086| . 

R. D. Young, D. B. Leinweber and A. W. Thomas, Prog. Part. Nucl. Phys. 50, 399 (2003) 
|arXiv:hep-lat /021203l1 . 

S. R. Beane, Nucl. Phys. B 695, 192 (2004) arXiv:hep-lat /0403030| . 

A. W. Thomas, P. A. M. Guichon, D. B. Leinweber and R. D. Young, Prog. Theor. Phys. 
Suppl. 156, 124 (2004) |arXiv:nucl-th/0411014| . 

D. B. Leinweber, A. W. Thomas and R. D. Young, Nucl. Phys. A 755, 59 (2005) 
|arXiv:hep-lat/05010 28|. 

J. F. Donoghue, B. R. Holstein and B. Borasoy, Phys. Rev. D 59, 036002 (1999) 
|arXiv:hep-ph/9804281| . 

D. Djukanovic, M. R. Schindler, J. Gegeha and S. Scherer, Phys. Rev. D 72, 045002 (2005). 
R. D. Young, D. B. Leinweber and A. W. Thomas, Phys. Rev. D 71, 014001 (2005) 
|arXiv:hep-la t /040 600T] . 

M. F. L. Golterman and K. C. L. Leung, Phys. Rev. D 57, 5703 (1998) ,arXiv:hep-lat/9711033l . 



S. R. Sharpe and N. Shoresh, Phys. Rev. D 64, 114510 (2001) arXiv:hep-lat /0108003| . 
J. W. Chen and M. J. Savage, Phys. Rev. D 65, 094001 (2002) |arXiv:hep-lat/0111050| . 
S. R. Beane and M. J. Savage, Nucl. Phys. A 709, 319 (2002) |arXiv:hep-lat /0203003|. 



29 



D. B. Leinweber, Phys. Rev. D 69, 014005 (2004) |arXiv:hep-lat/Q211Q17| . 

D. Arndt and B. C. Tiburzi, Phys. Rev. D 68, 094501 (2003) |arXiv:hep-la t /0307003| . 



D. Arndt and C.-J. D. Lin, Phys. Rev. D 70, 014503 (2004) |arXiv:hep-lat /0403012 1 . 

J. Bijnens, N. Danielsson and T. A. Lahde, Phys. Rev. D 70, 111503 (2004) 

|arXiv:hep-lat/0406017| . 

W. Detmold and C.-J. D. Lin, Phys. Rev. D 71, 054510 (2005) |arXiv:hep-la t /0501007" . 

J. Bijnens, P. Gosdzinsky and P. Talavera, Nucl. Phys. B 501, 495 (1997) 

|arXiv:hep-ph/9704212 . 

P. C. Bruns and U. G. Meissner, Eur. Phys. J. C 40, 97 (2005) |arXiv:hep-ph/0411223| . 



C. R. Allton, W. Armour, D. B. Leinweber, A. W. Thomas and R. D. Young, Phys. Lett. B 
628, 125 (2005) arXiv:hep-lat/0504022]. 

D. B. Leinweber and T. D. Cohen, Phys. Rev. D 49 (1994) 3512 |arXiv:he p-ph/930726T1 . 

D. B. Leinweber, A. W. Thomas, K. Tsushima and S. V. Wright, Phys. Rev. D 64, 094502 

(2001) :arXiv:hep-lat/0104013| . 

J. N. Labrenz and S. R. Sharpe, Phys. Rev. D 54, 4595 (1996) |arXiv:hep-lat /9605034| . 

C. K. Chow and S. J. Rey, Nucl. Phys. B 528, 303 (1998) |arXiv:hep-ph /9708432 . 

R. D. Young, D. B. Leinweber, A. W. Thomas and S. V. Wright, Phys. Rev. D 66, 094507 

(2002) !arXiv:hep-lat/02Q5017| . 



R. Sommer, Nucl. Phys. B 411, 839 (1994) ■arX iv:hep-lat/9310022| . R. G. Edwards, 
U. M. Heller and T. R. Klassen, Nucl. Phys. B 517, 377 (1998) |arXiv:hep-lat/9711003|. 



R. G. Edwards, U. M. Heller and T. R. Klassen, Phys. Rev. Lett. 80, 3448 (1998) 
|arXiv:hep-lat/9711052i . J. M. Zanotti, B. Lasscock, D. B. Leinweber and A. G. Williams, 
Phys. Rev. D 71, 034510 (2005) arXiv:hep-lat/0405015] . 
O. Bar, Nucl. Phys. Proc. Suppl. 140, 106 (2005) ■ arXiv:hep-lat/0409123| . 

G. Rupak and N. Shoresh, Phys. Rev. D 66, 054503 (2002) |arXiv:hep-lat /0201019| . 

O. Bar, G. Rupak and N. Shoresh, Phys. Rev. D 70, 034508 (2004) |arXiv:hep-la t /030602r . 
S. Aoki, Phys. Rev. D 68, 054508 (2003) |arXiv:hep-lat /030602f] . 

S. R. Beane and M. J. Savage, Phys. Rev. D 68, 114502 (2003) |arXiv:hep-lat /0306036| . 

H. R. Grigoryan and A. W. Thomas, |arXiv:hep-lat/05Q7028| 



The Adelaide approach also involves the A parameter which is discussed shortly. 



30 



